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Young star cluster evolution and metallicity 
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Abstract. Young star clusters (SCs) are the cradle of stars and the site of important dy- 
namical processes. We present W-body simulations of young SCs including recipes for 
metal-dependent stellar evolution and mass loss by stellar winds. We show that metallic- 
ity affects significantly the collapse and post-core collapse phase, provided that the core 
collapse timescale is of the same order of magnitude as the lifetime of massive stars. In 
particular, the reversal of core collapse is faster for metal-rich SCs, where stellar winds are 
stronger. As a consequence, the half-mass radius of metal-poor SCs expands more than that 
of metal-rich SCs. 
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1. Introduction 

The evolution of a star cluster (SC) is con- 
nected with the metallicity of its stars, be- 
cause the metallicity influences stellar winds 
(e.g. Kudritzki et al. 1987; Vink et al. 2001), 
remnant formation (e.g. Mapelli et al. 2009) 
and other properties of stars (e.g. Hurley et al. 
2000). Young (< 100 Myr) massive (> 10 3 
M ) SCs may have a particularly short two- 
body relaxation time 10 - 100 Myr, e.g. 
Portegies Zwart 2004). This implies that the 
core collapse and the post-core collapse phase 
occur on a timescale similar to the lifetime 
of massive stars. Thus, the peak of mass loss 
by stellar winds coincides approximately with 
the epoch of SC core collapse. The removal of 
mass by stellar winds and core-collapse super- 
novae (SNe) makes the SC potential well shal- 
lower, contributing to reverse the core collapse. 



In this paper, we discuss the results of V-body 
simulations of young SCs, including recipes 
for metal-dependent stellar evolution and stel- 
lar winds. We show that stellar metallicity can 
significantly affect the structural properties of 
SCs in the early post-core collapse phase. 

2. Simulations 

We perform V-body simulations of SCs us- 
ing the Starlab public software environment 
(Portegies Zwart et al. 2001). We modified 
Starlab, to include metal-dependent stellar evo- 
lution and recipes for stellar winds by Vink 
et al. (2001; see also Mapelli et al. 2013 for 
more details on the code). We simulated young 
intermediate-mass SCs, generated according 
to a multi-mass King model, with total mass 
M T0 T = 3000 - 4000 M , initial core radius 
r c = 0.4 pc, concentration c = 1.03. The stars 
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Fig. 1. Top: core radius r c as a function of time t. 
Bottom: half-mass radius rh m as a function of t. Red 
solid line: metallicity Z = 0.01 Z Q ; black dashed 
line: 0.1 Z Q ; blue dotted line: 1 Z Q . Magenta line in 
the bottom panel: analytical prediction (r hm oc t 2/3 ). 

in the SC follow a Kroupa (2001) initial mass 
function. We consider three different metallic - 
ities: Z = 0.01, 0.1 and 1 Z . We ran 100 real- 
izations of the same SC (by changing the ran- 
dom seed) for each metallicity, to filter out the 
statistical fluctuations. 

3. Results and conclusions 

The core collapse occurs at t ~ 3 Myr in all 
simulations (top panel of Fig. 1). The reversal 
of core collapse depends strongly on the stel- 
lar metallicity, as the post-collapse r c is a fac- 
tor of 1.5 larger in metal-rich SCs than in 
metal-poor SCs. The reason is that mass loss 
by stellar winds and SNe is stronger in metal- 
rich SCs, making the core potential well shal- 
lower. At later times, the core radii at differ- 
ent Z become similar again, as the effect of 
stellar winds is over, and the evolution of the 
core is completely determined by three-body 
encounters. The half-mass radius (rh m , bottom 
panel of Fig. 1) starts expanding after core col- 
lapse, according to the well-known analytical 
model (rh m 00 t 2/i , e.g. Elson et al. 1987). On 



the long-term evolution (f ~ 100 Myr), the 
half-mass radius in metal poor SCs is ^ 10 
per cent larger than in metal-rich SCs. The 
reason is that, in metal-poor SCs, the rever- 
sal of core collapse is slower, implying higher 
core densities for a longer time. This means 
that the rate of dynamical interactions in the 
core of metal-poor SCs is higher, and more ki- 
netic energy is pumped into the halo, increas- 
ing rhm. This result is in agreement with the 
recent simulations by Schulman et al. (2012). 
Other studies (Downing 2012; Sippel et al. 
2012) do not find important differences in rh m , 
as they consider systems with much longer re- 
laxation time. Our results open interesting per- 
spectives on the study of the dynamical evolu- 
tion of young massive SCs. 
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